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Abstract 



We present a complete calculation of the K13 decays K'^ n^l^vg and n^t^vg to 0{p^, [nid — rau)p'^ ^e^p^) 

in chiral perturbation theory with virtual photons and leptons. We introduce the concept of generalized form factors 
and kinematical densities in the presence of electromagnetism, and propose a possible treatment of the real photon 
emission in K^^ decays. We illustrate our results by applying them to the extraction of the Kobayashi-Maskawa matrix 
element |Kts| from the experimental K'^^ decay parameters. 
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1 Introduction 

Semileptonic kaon decays have played a central role in 
our understanding of flavour physics. In particular, the 
JTeS decay mode is often presented as the most accurate 
and theoretically cleanest source for the extraction of the 
Kobayashi-Maskawa matrix element 0- However, 

as new high precision experiments are planned and the 
theoretical tool of effective field theory has been pushed to 
higher orders, an update of the theoretical analysis of such 
decays is needed. In particular, it is interesting to improve 
as much as possible the theoretical analysis underlying the 
extraction of \Vus\ - This is presently based on the following 
ingredients: a calculation of the hadronic form factors 
at order in chiral perturbation theory (CHPT) 
supplemented by a model-dependent estimate of the order 

effects As for the radiative corrections, the model 
independent short distance leading logarithms [^|^ have 
been included, together with a model dependent estimate 
of the long distance contributions P,p|,p^,pl|. 

It is the purpose of this paper to work out the radiative 
corrections to the if ^3 decays {£ — e, ^) within the frame- 
work of CHPT, the effective theory of the standard model 
at low energy. The appropriate formalism for including vir- 
tual photons in purely mesonic low-energy processes has 
been presented in This scheme was then ex- 

tended for the treatment of semileptonic interactions by 
the additional inclusion of virtual leptons The goal 
of our analysis is to obtain the size of such corrections 
and quantify the theoretical uncertainty to be assigned 
to them. This will allow us to understand how well one 
can know \Vus\ once new high statistics experiments col- 
lect data and the full order CHPT analysis becomes 
available JlGf . 

The outline of the paper is as follows. In Sect. ^ we 
recapitulate CHPT in the presence of virtual photons and 
leptons. The basic Ki^ phenomenology in the absence of 
electromagnetism is reviewed in Sect. 0. The structure of 
the radiative corrections to Ki^ decays is discussed in Sect. 
^. The Ki3 amphtudes to O(p^) and C'(e^p^) are calcu- 
lated in Sect. |[ In Sect. ^ we propose a specific treatment 
of the real photon emission in the K'^^ case. In Sect. ^ 
we illustrate our general considerations by a numerical 
study of the K^^ decay and the description of a procedure 
to extract the Kobayashi-Maskawa matrix element \Vus\ 
from experimental data. Our conclusions are summarized 
in Sect. ^ Appendix ^ contains a summary of integrals 
appearing in various mesonic one-loop amplitudes and Ap- 
pendix B collects several photonic loop functions. In Ap- 
pendix C we report a set of coefficients appearing in the 
next-to-leading order expansion of the form factors. 



2 The standard model at low energies 

The appropriate theoretical framework for the analysis of 
electromagnetic effects in semileptonic kaon decays is a 
low-energy effective quantum field theory where the asymp- 
totic states consist of the pseudoscalar octet, the photon 



and the light leptons . The corresponding lowest-order 
effective Lagrangian is given by 



off 



f {u^u^^ + x+) + e^F^ZiQ^QT) - -.F^uF^" 



e 



[£i^ 



(2.1) 



The symbol ( ) denotes the trace in three-dimensional 
flavour space, and 



(2.2) 



The photon field and the leptons £, vg {£ — e, /i) are 
contained in (2.2) by adding appropriate terms to the 
usual external vector and axial- vector sources u^, a^: 



(2.3) 



The 3x3 matrices Qf\, Ql spurion fields. At the 
end, one identifies Ql"r with the quark charge matrix 



2/3 
-1/3 
-1/3 



g"" = 

whereas the weak spurion is taken at 
Ql = -2V2 Gf 



Vud Vus 






(2.4) 



(2.5) 



where Gp is the Fermi coupling constant and Vud, Vus are 
Kobayashi-Maskawa matrix elements. For the construc- 
tion of the effective Lagrangian it is also convenient to 
define 



(2.6) 



_Fo denotes the pion decay constant in the chiral limit and 
in the absence of electroweak interactions. Explicit chiral 



symmetry breaking is included in x+ = m + 
where x is proportional to the quark mass matrix: 



jUu 

TUd 

m, 



ux^u 



(2.7) 



The factor Bq is related to the quark condensate in the chi- 
ral limit by (0|^q|0) = —FqBq. This leads to the (lowest- 
order) expressions for the pseudoscalar masses 



2Bom + 2e^ZFl, 



= 2Bom, 



M 



Bo 



{rris + fh) 



2£(2) 



{tTLs — to) 



=Bn 



"-(-) 
K 



V3 

(TOs + to) H -^(ms — m) 



V3 



4 / TO 

3B0 (to. + - 



+ 2e^ZF^, 



(2. 



3 



The tree-level mixing angle e^^-' is given by 
V3 ma-mu 



M - 



(2.9) 



the symbol fh stands for the mean value of the light quark 
masses, 

m = + TOd). (2-10) 

For later use, we also denote the isospin limits (m„ — nid, 
e = 0) of the pion mass and the kaon mass, respectively, 

by 

= 2Bam, M|. = Ba{m, + m). (2.11) 



Using (2^), the numerical value of the coupling constant 
Z ~ 0.8 can be determined from the mass difference of 
the charged pions. 

In our calculation we include terms in the low-energy 
expansion up to order p'*, (m„ — md)p'^ and e^p^. 

The most general local action at next-to-leading order 
can be written as the sum of four terms, Cpi + LfS^i -I- 



lept 



L^. The first one, H^i is the well-known Gasser- 



Leutwyler Lagrangian [| in the presence of the generalized 
external sources introduced in (2.3): 



£„4 = 



+ L3 (Uf,u>'u^u'') + Li {u^.u^') (x+) 

+ i(2L8 + Li2)(x+) + J(2L8 - ii2)(X-) 



with 



^(iio — 2Lll)(/-p,y/ 



(2.12) 



(2.13) 



The second term, Ce^p2 , denotes the interaction of dynam- 
ical photons with hadronic degrees of freedom |2[|l|,|l|. 
Its expression is: 



£g2p2 



^0 ((sr 



- K3 [{Qi:-u, 

+ K, ((srsr + QTQrw^n 
+ mrf + {QT?) (x+) 

+ ^9([(Qr)' + (Qr)']x+) 

+ ^10 {{QrQT + QTQr)x+) 



i^ii ((srsr - srsDx-) 
iKi2 ([(v^QDQr - srv^Qr 

i^i3 ((v^Qr)(v^Qr)> 
i^i4 ((v^Qr)(v^Qr) 

+ (v^Qr)(v'^sr)>V (2.14) 



where 



with 



(2.15) 



(2.16) 



The presence of virtual leptons requires also the addition 
of the "leptonic" term (l5| 



^Icpt — 



{Fo^[x,h,,yeL{unQ'^^,QU) 

+ iX5h^,,eL{Q'[^''QT) + h.c. 
+ Xeiii ^ + e4)£ 
+ X7meUy 



(2.17) 



In >Cicpt we consider only terms quadratic in the lepton 
fields and at most linear in Gp. Finally, also a photon 
Lagrangian 



= e^XsF^.F^'', 



F^u ^ d^A, - d,A^, (2.18) 



has to be added |15 



The low-energy couplings Li, Ki, Xi arising here are 
divergent (except L3, L7, K-^, K13, Kn and Xi). They 
absorb the poles of the one-loop graphs via the renormal- 
ization 



= L[(Ai) + r,yl(Ai), i = l,...,12, 
K,^Kl-{fi) + S,A{fi), i = l,...,14, 
X,=X:ifi)+SM, z = l,...,8. 



(2.19) 



with 

1(f) 



20) 
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in the dimensional regularization scheme. The coefRcients 
Fi and Si can be found in [Q and in respectively. 
Their values are not modified by the presence of virtual 
leptons as long as contributions of 0{Gp) are neglected. 
The "leptonic" coefficients Si have been determined rather 
recently in by using super-heat-kernel techniques . 

In order to match our low-energy effective theory to 
the standard model of strong and electroweak interactions, 
we have to specify the precise physical meaning of the 
parameter Gp. In the presence of electromagnetism, this 
identification is somewhat ambiguous. The shift 



Gf -> Gf(1 + e^S) 
induces the change (see Eq. (4.8) of p5[) 



p2 
-^0 



i^) -2e2(5^^(^^- 



corresponding to 



Xq — > Xq — 25, Xj — > Xy + 26. 



(2.21) 

mi)e, 
(2.22) 

(2.23) 



3.1 Invariant amplitude 

We shall consider the generic Ki^ process 

K{pk) 7r(p7r) t'^ivi) vi>{pu)- 
The invariant amplitude reads 



where 

;^ = S(P.)7^(1-75)«(P£), c 



(3.1) 



(3.2) 



1 for 
— for 



The expression in parentheses corresponds to the matrix 
element {n{pTr)\V^~^^\K{pK)), expressed in terms of the 
form factors f^"^ (t) . The hadronic form factors depend on 
the single variable t — {pK — p-kY . 



In other words, some part of the electromagnetic contri- 
butions may always be shuffled from Xq to Gp or vice 
versa. (The coupling constant X-; does not appear in ob- 
servable quantities as it is always absorbed by the mass 
renormalization of the charged leptons.) 

Following we identify Gp with the muon decay 
constant. To order a, Gp can be related to the measured 
muon decay width by |lq] 



r{ii^ aU) 



1927r3 



(2.24) 
With this choice 



with /(x) = 1 - 8a; - Ux^ In x 
of Gp, the (universal) short-distance electromagnetic cor- 
rection 1^,01 of semileptonic charged current amplitudes 
is fully contained in the coupling constant Xg. To display 
this dependence explicitly, we pull out the short-distance 



part Xq^ by the decomposition 



xii^,) = xt"" + xl{^il 



(2.25) 



where 



^ = l-S^ 



EW 



{Mp,Mz), (2.26) 



which defines |j] also the short-distance enhancement fac- 
tor S'ESN{Mp,Mz) (to leading order). With this splitting 
of Xg, we expect its "long-distance part" Xg (Mp) having 
the typical size of a low-energy constant in CHPT. 



3.2 Dalitz Plot density 

It is customary to analyze the spin-averaged decay dis- 
tribution p{y,z) for K^'i. It depends on two variables, for 
which we choose: 



2p7r • PK 2Et, 



Ml 



M, 



K 



y 



2pK ■ Pi 

Mi 



where E^^ {Ei) is the pion (charged lepton) energy in the 
kaon rest frame, and Mk indicates the mass of the de- 
caying kaon. Alternatively one may also use two of the 
Lorentz invariants 

t = (pk -Pt,)"^, u={pK-pef, s = (p^ +pt)'^. {3A) 

Then the distribution (without radiative corrections) reads 

(3.5) 

with 



^=^^ 128.3 
Moreover, defining 



,2G||K.pAf| 



J dydz p^''Hy,z). (3.6) 
_rnl _Ml 



Ml 



Ml' 



3 Basic Kis phenomenology 

In this section we report the basic formulae of K13 phe- 
nomenology in the absence of radiative corrections |^. 

^ In the present section Mk and stand for the physical 
masses of the corresponding mesons involved in the process. 



the kinematical densities are: 

Af\y, z) = A{z + y- 1)(1 - y) + r,{Ay + 3z - 3) 
- 4r^ + ri{r^ - ri), 



4°^(y,z) = 2re{3-2y-z + n~r^), 
A'f \y,z) = reil + r^, - z ~ re). 



(3.8) 
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In the analysis of K^^ decays, the terms proportional to 



Aj'^g can be neglected, being proportional to 
The physical domain T) is defined by 

< 2/ < 1 + r£ - r^, 
a(2/) - Ky) <z< a{y) + 



-6 



10 



(3.9) 



where 

b{y)-- 

or, equivalently, 



(2 -y)il + n + r^- y) 
2(1 + r, - y) 



^j/^ - 4r£ (1 + - - j/) 
2{l + rt~y) 



(3.10) 



2\/?V < 2 < 1 + Ttt - Tf, 

c(z) - d{z) <y< c{z) + d{z), 



(3.11) 



where 



c(z) 
d(z) 



(2 - z) (1 + + r£ ^ z) 

2(l + r^-z) 
\/z2~^4rv (1 + r^r - r£ - z) 



2(1 



(3.12) 



4 Structure of radiative corrections to 
decays 

In decays involving charged particles, the observable quan- 
tity always involves an inclusive sum over the parent mode 
and final states with additional photons. Therefore, when 
considering radiative corrections one must include the ef- 
fect of real photon radiation as well as virtual electromag- 
netic corrections. Moreover, from a theoretical point of 
view, only such an inclusive sum is free of infrared (IR) 
divergences order by order in a. In the present section we 
give an overview on the analysis of Ki^ decays with in- 
clusion of radiative corrections. We show how electromag- 
netic corrections can be accounted for by using generalized 
form factors and kinematical densities, whose specific form 
we describe in detail in the following sections. 

The virtual corrections will in general produce a shift 
in the invariant amplitude A4 ^ Ad + AM. . On the other 
hand, the radiation of real photons (governed by an am- 
plitude M'') will produce a shift to the decay distribution 
and rate. 

The virtual electromagnetic correctio ns do not alter 
the structure of the invariant amplitude ( |3.2| ) in terms of 
the form factors f± (t) , but change the form factors them- 
selves. We denote by F±{t, v) the full form factors includ- 
ing virtual electromagnetic corrections. Note that such ob- 
jects now depend also on a second kinematical variable 
(this is because F± cannot be interpreted anymore as ma- 
trix elements of a quark current between hadronic states). 
The variable v is taken as u = {pk — Pif' for K^^ and 
s = (ptt + Pe)'^ for X^g. We observe here that F± contain 



infrared singularities, due to low-momentum virtual pho- 
tons: we regularize them by introducing a small photon 
mass Mj. Summarizing, the effect o f vi rtual corrections is 
to change the invariant amplitude (3^) as follows: 



M [ft'\t)J^\t)]-^M[F+it,v),F^it,v)]. (4.1) 

Moreover, it is convenient to factor out of F± the long dis- 
tance component rc{v , rnj , A-P ; M^) of loop amplitudes, 
which generates infrared and Coulomb singularities, as fol- 
lows: 

F±it,v)=\l + ^r,{v,m'i,M^;M^)\ U{t,v). (4.2) 

Fc depends on the mass mg of the charged lepton, the mass 
M of the charged meson, the variable v and has a logarith- 
mic dependence with respect to the infrared regulator M^. 
The definition of F^ is not unique (due to possible shifts 
i n th e infrared finite parts) , and we give o ur d eifinition in 
(5J_). Once one specifies the form of I^;, (O) serves also 
as a definition of the structure-dependent effective form 
factors f±{t,v). In the next section we shall present the 
calculation of f± at next to leading order in CHPT. 

Let us now comment on the radiative amplitude M'^ , 
and how it affects the decay distribution. In the intermedi- 
ate stages of this work we shall use Low's theorem jl^ to 
obtain the leading components of M'^ at small photon mo- 
mentum. This is consistent with a full calculation at order 
^2p2 CHPTQ. However, it allows us to be more general 
in the derivation of the radiatively-corrected decay dis- 
tribution, by including at intermediate steps some higher 
order terms in the chiral counting. We restore the correct 
chiral counting by inserting in the final expression for the 
decay rate the form factors as calculated at 0{p'^, e^p^) in 
CHPT. Let us now sketch the analysis. The radiative am- 
plitude can be expanded in powers of the photon energy 



M'' 



where 



-■^(0) 



El 



(4.3) 



Gauge invariance relates Mj_.y^ ■"'i(o) 



(4.4) 

and ^^(0) ^^'^ non- 
radiative amplitude Al, and thus to the full form factors 
^±(^1^^)- Upon taking the square modulus and summing 
over spins, the radiative amplitude generates a correction 
p^(y, z) to the Dalitz plot density of ( [3.5[ ). The observable 
distribution is now the sum 



p{y,z) ^ p^"^y,z) + p^{y,z). 



(4.5) 



Both terms on the right hand side of this equation depend 
on the full form factors F± and contain infrared diver- 
gences (from virtual or real soft photons). Upon combin- 
ing them, the observable density can be written in terms 
of new kinematical densities Ai and the effective form fac- 
tors f±{t,v) defined in (jij), 

p{y,z)=Af Sew{Mp,Mz) (4.6) 

X 'A^\f+{t,v)\' + A2f+{t,v)f-{t,v) + A^\f^{t,v)\'^ 



^ At this order CHPT reproduces the results of Low's theo- 
rem, with f+{t) = 1 and f-{t) — 0. 
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where we have pulled out the short-distance enhancement 
factor S'ew discussed in Sect. ^j. To first order in a, the 
kinematical densities are given by 

My,z) = Af\y,z) [l + A'^{y,z)\ + A'^{v,z). (4.7) 

The function Z\^^(j/, z) arises by combining the contribu- 
tions from 

? and rc{v, ■mj,M'^; M^). Although the 
individual contributions contain infrared divergences, the 
sum IS finite. The factors Z\P(?/, z) originate from averag- 
ing the remaining terms of |A^'''p [see (i.3)] and are IR 
finite. Note that both A^^{y, z) and /\-^, z) are sensi- 
tive to the treatment of the real photon emission. Details 
on these corrections are given in Sect. |[ 
Let us finally note that, in principle, the radiative ampli- 
tude generates new terms in the density, proportional to 
derivatives of form factors. These terms would only arise 
at order e^p* and higher in CHPT, and therefore we have 
suppressed them in (4.6). 

Summarizing, the resulting Dalitz density is formally 
identical to the unperturbed one, but now one has 



Af\y,z)^My,z). 



(4. 



We want to stress here that although the structure de- 
pendent electromagnetic effects in f±{t,v) are due to the 
interplay between QCD dynamics and QED, the correc- 
tions of order a included in the densities Ai{y,z) are 
universal in the sense that they are fixed by gauge in- 
variance and kinematics. However they do depend on the 
choice of the experimental cuts to the photon spectrum 
and they exhibit a certain ambiguity due to possible dif- 
ferent defi niti ons of the function rc{v, mj, M^; M^) intro- 



duced in (4.2) (we adopt the form given in ( p.l| ) below). 
We suggest that these modified densities be used in the 
data analysis. This does not introduce any model depen- 
dence, and takes care of very long distance electromagnetic 
corrections. Such an experimental analysis of the Dalitz 
plot density could then provide valuable information on 
f±{t,v), to be confronted with theoretical calculations of 
pure QCD and QCD+QED effects. 



5 K^s amplitudes at order 
CHPT 



and 



In this section we analyze the form factors relevant for 
decays in the framework of CHPT, including elec- 
tromagnetic corrections at order e^p^. As for the non- 
electromagnetic part, we shall give the results at order 
in standard CHPT, including the effect of strong isospin 
breaking (to„ rad). In principle, the order corrections 
can be easily included in our formalism, in order to make 
a more accurate phenomenological analysis. 

In Sect. 1^ we have introduced the effective form factors 
f±(t,v), from which the long distance virtual photon ef- 
fects have been removed [see (i^)]. We also stressed that 
in order to give a meaningful definition of such objects one 



has to specify the function rc{v , rnj , M"^ ; Mj) , arising in 
perturbation theory from photonic loop contributions to 
the form factor Fj^.{t,v): 



r,{v,mlM^;M^) = 
2M'YC{v,mlM') + 2log^(^l 



XYlogX 



(5.1) 



In Appendix |b| we define the auxiliary variables X, R, Y, 
as well as the function C{v, mjjM"^). encodes the lead- 
ing Coulomb interaction between the charged particles in- 
volved in the decay {K^ and for if^, 7r~ and for 
K^^). It is IR divergent and singular when the two charged 
particles are relatively at rest. In the case of if + decays, 
this singularity is outside the physical region, while it oc- 
curs on its boundary for the K'^ decay, implying a larger 
overall effect. It is convenient to write the form factors in a 
form suitable for a resummation of (potentially large) IR 
and Coulomb effects. At the same time we want this rep- 
resentation to reproduce the correct CHPT result at order 
e^p^ . We achieve this by writing the full form factors as 
follows, in terms of the effective form factors: 



F±{t,v) = [1 + —r,{v,ri4,M^-M^)\ f±{t,v) , 



Let us now comment on each term appearing in ( p.2[ ): 

^ /±(i) represent the pure QCD contributions to the 
form factors (in principle at any order in the chiral ex- 
pansion), plus the electromagnetic contributions up to 
order e^p^ gen erated by the non-derivative Lagrangian 
defined in (2T) (proportional to Z). Diagrammatically 
they arise from purely mesonic graphs. In the defini- 
tion of /f'^ (t) , we have included also the electromag- 
netic counterterms relevant to ir^-rj mixing [see (^.4|) 



and (5.5) below] 



jEM loc represent the local effects of virtual photon 
exchange. 

yEM loop ^^^^ represent the non-local photonic loop con- 
tribution (once the Coulomb term has been removed) . 

We have checked that in all relevant cases ^°°^{v) 
are smooth functions of the variable v in the physical 
region, allowing one to perform a linear expansion in 
the Dalitz variables y, z. 



Each term in the above decomposition is scale- independent. 
We now give the full expressions for the form factor cornp- 
ponents, in terms of loop functions defined in Appendix H. 
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5.1 The form factors f^^'^°{t,u) 



For the mesonic contribution is given by: 

1 3 



1287r2F2 



M2 

(3 + V3s(2))M,^ln^ 



2(3- V3e(2))M|.o In 



V3e(2) 



(5.3) 



This expression is essentiahy the pure QCD form factor 
at 0{p'^) [||, except for the inclusion of electromagnetic 
contributions to the meson masses and for the additional 
contribution £^^1' ^? Tr^-ry mixing at ©(e^p^J^ jlJ 



The function HpQ(t) is reported in Appendix and 
the leadi ng o rder ir^-rj mixing angle e'-^^ has alr eady been 
given in (2.9). The remaining quantities in (5.3) are given 

by 



-(4) 



2e(2) 



3(4^Fo)2(M2-M2) 
X <i {Airf 64 [3^7 + Ll{^x)] [Ml ~ M^f 



M. 



Ml 



- M^iM], - Ml) log + Ml{Ml - iMl) log 

- 2Ml{Ml - 2Ml) log ^ - 2MUmI - M^^ 

(5.4) 

and ^ 

(4) _ 2v/3aM|- 
^EM- i08 7r(Af2-M2) 

X (2(4^)2 [ - QKl{^l) + 3A-I(/i) + 2i^5^(A*) + 2Kl{^l) 



-9Z(log^ + l 



They are related to the one-loop off-diagonal element of 
the squared mass matrix in the Tr^-yy sector [031 by 



-(4) , ^(4) 



EM 



(M2-M2)- 



(5.6) 



The analogous contribution for the /_ form factor is given 
by: 

/_(t) = ^(l + ^)(Mi-M2) 



V3' 



2567r2 



In- 



M' 



3 The last -1 in Eq. (5.21) of should be replaced by -1-1; 
the preceding formula (5.20) is correct. 



M2 

Mi: 

2(3-N/3e(2))M2± In^ 
Ml. 



+ (l + 3V3e(2))Af2„ in—i 



A 



PQ, 



PQ 



2t 



bpQit) KpQ{t) 



(5.7) 



The sum in the last line runs over meson pairs occurrinj 
in loop diagrams. The functions KpQ{t) and Hpgit) ^ 
are reported in Appendix while the relevant coefficients 
apqit) and hpQ{t) are listed in Appendix |c[ 

The local electromagnetic contributions are given by: 



/EM-loc ^ 



327r2 



3 -I- In 



31n- 



M 



(5.8) 



/EM-loc ^ 



2/^3^ (^i) - Kl{p.) ~ -{K^if,) + Klip)) 



+ X, - xUf^) + x^ifi) 



327r2 V 



1 - 21n 



(3+—) In 



(5.9) 



Finally, the loop contributions are given by 



j,EM-loop. N 



47r 



ri{u,mlMl) ± r2(u,m^,M|.)J , 

(5.10) 



(5 5) with ri^2iv, mj, M ) given in Appendix | 



5.2 The form factors /^°^ (t, s) 

In the case of f^^^^ , the mesonic loop contribution is 
given by 

~ 1 3 

/+(*) = 1 + 2^K+TT'>{t) + -^Hj^+jjit) + HkOtt- (t) 



+ ^e^^^ [HK.it) -H Knit)] , 



(5.11) 



and 



4 2e(2) 

l(l + i^)(M^-M^2 

■3 M2 



i5(M) 



2567r^ 



In- 



8 



1287r2i^2 



(3 + 2V3£(2))M2ln. " 



M ± 



^" - (3 + 2V3£(2))M^oln 



PQ 



2t 



dpqit) KpQ{t) 



+ dpQ{t)^HpQ} 



(5.12) 



As previously, the sum runs over meson pairs occurring 
in loop diagrams. The coefficients cpQ{t) and dpQ{t\ for 
the relevant meson pairs, are reported in Appendix O. 
The local electromagnetic terms are given by: 



/■EM-loc 



Ana 
1 



327r2 



In 



M2^ 



3 In ■ 



,(5.13) 



6 Corrections to kinematical densities 

As an illustration of the general considerations of Sect. 
^, we present here in detail a possible treatment of the 
contribution of the real photon emission process 

K+{pk) ^ Ti\p.)l+{pi)Mpu)-i{p^). (6.1) 

Following the procedure proposed by Ginsberg |^ , we de- 
fine the kinematical variable 



X= {p^+ p~ff = {pK - Pit- Ptf 



(6.2) 



A possible choice of the purely radiative part of the de- 
cay distribution (corresponding to a specific analysis of 
the experimental data) is to accept all pion and charged 
lepton energies within the whole K^s Dalitz plot T) and 
all kinematic ally allowed values of the Lorentz invariant 
X defined in (x2). (The variable x determines the angle 
between the pion and lepton momentum for given energies 
£'£.) This translates into: 



jEM-loc ^ 



1 - 2 In ■ 



327r2 



(3 



Z ^ 



In ■ 



M2 



(5.14) 



with 



pol 



(6.3) 



Finally, the loop contributions are: 



f± ^(s) = ^ [r2{s,mlMl) ± rr{s,mlMl) 



(5.15) 



5.3 The effective form factors for Kez decays 

In the analysis of K^-i decays (where only _F1|_ is observ- 
able), in order to avoid unnecessary complications we can 
write 



U{t)[ 



a 



l + __r,(.,^2^M2;M,)+/™-P 



with 



EM-loc 



(«)}, 
(5.16) 

(5.17) 



Exploiting the ambiguity in the definition of the function 
Fc , here we choose to shift a smooth function of v from the 
form factor to Fc itself. The advantage is that the effective 
form factor becomes a function of the single variable t. 
This trick is not possible in the treatment of -fT^a decays, 
as the loop contributions f+{v) and /-(w) are different. 



1 r 



yz + ^(j/2 - Are){z'^ - 4r^) 



(6.4) 



In (6^) we have extended the integration over the whole 
range of the invariant mass of t he unobserved Vi 7 system. 
The integrals occurring in (|6.3D have the general form M 



Im,niPl,P2; P, My) 

1 r d^q d^k 
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5(4) (p. 



k) 



fco (pi • fc + M2/2)"(p2 • fc + M2/2)' 



(6.5) 

The results for these integrals in the limit = can be 
found in the Appendix of ( We h ave checked these ex- 
pressions.) U sing the definition (6^), the radiative decay 
distribution (6.3) can be written as M 



K 



1287r3 



j dx'^Cm,nIm,n{Pl,PK]PK -P-^ ~Pi,0) 



(6.6) 



9 



M2 



where the infrared divergences are now confined to 
Io{y,z;M^) = 



dx - 2pK ■ pih^iipi, -pk;Pk ~Ptt - Pe; 

- mll2fi{pt,-pK\PK ~ Pi, ~ Pf,M^) . (6.7) 



The explicit form of the function Iq can be found in Eq. 
(27) of 1^. The coefficients c,n,n were given in Eq. (19) of 
P]. Note however the misprint for the values of c_i^o a-nd 
C1.-2 (see Erratum of 1^). (We have also checked this list 
of coefficients.) 



The function introduced in (4.7) can now be re- 
lated to Iq by 



Z\i«(y,z) 
where 



a 



Io{y,z;M^ 



1 



reiu,mlMl;M^) 



(6.8) 



(6.9) 



j.,^/r. + ij/™-P for^ = e 

for ^ = ^ 



An analytic expression of the integral occurring in ( |6.6[ ) 
was given in Appendix B of JTTI in terms of the quantities 



dx c 



2^0 



(6.10) 



Note that the quantity JQ[i) given in Eq. (A9) of |^ 
(which is needed for the evaluation of U^) contains two 
crucial mistakes: the plus-sign in the last line of (A9) 
should be replaced by a minus-sign, and 1/3,™'^'^! at the 
end of the first line of (A9) should simply read (3f^'^^. The 
second error is irrelevant for H — e but has disastrous con- 
sequences in a certain part of the Dalitz plot for ^ = /i. To 
the best of our knowledge none of these errors has been 



reported in an Erratum of 1 1 1 



The functions Z\J^ introduced in (4/7) can now be writ- 
ten as 



IB 



A- 



^2 



2a 



7 



2a 



7 



(6.11) 



where the symbol ^ used in |Tl[ stands for the ratio /_ / /_| 



7 Application to K^^ decay 

In this section we provide an illustrative analysis of the 
K^^ decay with inclusion of isospin-breaking and radiative 



corrections. The aim here is to illustrate how one might 
proceed in order to extract |V^s| rather than to give final 
numbers, which would anyhow still depend on the size of 
the hitherto unknown two-loop strong interaction contri- 
butions. Once they become known, these corrections could 
then easily be incorporated into the analysis presented 
here. 

As usual we provide a linear expansion for the effective 
form factor ^ (t), which reads: 



with 
and 



/r""(o)-/+(o)+/f 



df+it) 



EM-loc 



dt 



(7.1) 

(7.2) 
(7.3) 



Using the linear expansion for the form factor (t), 
the decay rate 



r{K+ 



7r"e+^.,(7)) := 

can be expressed as 

r{K+ ^ 7r"e+z.e(7)) = 



^"e+i/eT) (7.4) 



where 



N S^^iMp.Mz) |/f^"°(0)|2 /(A+), 



/(A+) = / dydzAi{y,z) 



(7.5) 



= flO + «1 + «2 ^^4 



(7.6) 



In ( |7.5D we have pulled out the short distance enhance- 
ment factor S-e\n (Mo, Mz) from the low energy constant 
X^ifi) [see dH)]. 

In order to extract \Vus\ we have to provide a theoret- 
ical estimate of the form factor at i = and the phase- 
space integral. 



7.1 Numerical estimate of / 



(0) 



In our numerical analysis, we use the physical meson masses 
0. We now describe the input used for the other quanti- 
ties occurring in the form factor ex pan sion. 

For the mixing angle defined in (2.9), we are using [pO[: 



£(2) = (1.061 ± 0.083) X 10" 



(7.7) 



The com bina tion of coupling constants 3L7 + L'^{Mp) 
entering in ( ^.4[) can be determined from two observables 
[Q: the deviation from the Gell-Mann--Okubo Mass for- 
mula. 



A 



GMO 



4Af|. - Af 2 - 3A/2 
A'P - Af2 ' 



(7.8) 
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which is well under control, and a quantity Am related to 
the 0{p'^) contributions for the ratio 



Mi 



nis + m 



{1 + A 



M2 2m 
Combining Eqs. (10.10) and (10.11) of 



M) 



(7.9) 



3LT + Ll{Mp 



A. 



GMO 



one finds 



24(M2 - M2) 32(Mi - M^) 
12(M2-Af2)2 /' 



where 



Using 



0.065 ±0.065 



from Leutwyler's analysis |22l, (7.10) implies 
3L7 + Ll{Mp) = (-0.33 ± 0.08) X 10~^ 



(7.10) 

(7.11) 
(7.12) 

(7.13) 



with a relatively small error. (Allowing also for higher 
order corrections one might prefer [|l| a more generous 
bound of, say, |Z\m| < 0.2 corresponding to SL^+LKAIp) — 
(-0.25 ±0.25) X 10-3.) 

The relevant combin atio n of electromagnetic coupling 
constants appearing in (5.5) has been estimated by Bij- 
nens and Prades ||22| . For our numerical analysis we add 
here an error of the typical size l/(47r)2 ~ 6.3 x 10"'^: 



X'-(Mp) := (-6X3 + + 2K5 
= (5.7 ±6.3) X 10"^ 



2K^Y{Mp) 



(7.14) 



In this way we find the following mesonic contributions to 
the form factor: 



/+(0) = 1.0002^1 + 0.0228- 



0.0055 • 



3L7 



- 1.06 X 10-2 
1.06 X 10-2 
hig(A/p) + 0.33 X 10-3 



0.0002 



K-{Mp] 



-0.33 X 10-3 
- 5.7 X 10-3 



5.7 X 10-3 
1.0002 ± 0.0018 ± 0.0013 ± 0.0002 
1.0002 ± 0.0022. 



(7.15) 



In the last line the three individual errors have been added 
in quadrature. The main contribu tion to the fina l error 
comes from the uncertainties in (7/7) and ( [7.13| ). The 
electrom agne tic contribution to 7r° — 77 mixing, entering 
through (5.5), has very little numerical impact on the cen- 
tral value and the uncertainty. 

For the coupling c ons tant 7^12 entering in the purely 
electromagnetic part (^) we use a value extracted from 



(7.16) 
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K{^{Mp) ^ (-4.0 ±0.5) X 10-3. 



Table 1. Coefficients entering the phase space integral 





ao 


0.1 


a-2 


a = 


0.09653 


0.3337 


0.4618 


a / 


0.09533 


0.3287 


0.4535 



For the (unknown) "leptonic" constants we resort to the 
usual bounds suggested by dimensional analysis: 

|Xe^(Afp)| < 6.3 X 10-3. (7.17) 

Eventually we find: 

Kl^{Mp)^Ax 10-3 



jEM-loc ^ 0.0032 - 0.0007 • 



0.0015- 



4 X 10-3 

XliMp) 



6.3 X 10-3 
= 0.0032 + 0.0016. 



0.0003 • 



6.3 X 10-3 

(7.18) 



In this contribution, the sizeable relative uncertainty is 
almost exclusively due to the poor present knowledge of 
Xi. Despite this, in the final result for /:f^^^°(0) this is an 
effect of 0.16%. 

Combining the results given above we obtain the final 
value: 

(7.19) 



f^^^'^iO) = 1.0034 ±0.0027. 



7.2 The phase space factor 

The theoretical prediction for the slope parameter is de- 
termined by the size of the low-energy constant Lg. With 



Ll{Mp) = (6.9 ±0.7) X 10-^ 



(7.20) 



we find 



0.0328 + 0.0321 
0.0328+0.0033. 



Ll{Mp) - 6.9 X 10- 
6.9 X 10-3 



(7.21) 



Here one can note the relatively large uncertainty induced 
by Lg(Mp). This suggests to use the measured value for 
the slope when evaluating the phase space integral. 

The numerical coe fficients ao,i.2 entering in the phase 
space expression (7.6) are reported in Table ^. The num- 
bers given here correspond to the specific prescription for 
the treatment of real photons described in the previous 
section: accept all pion and charged lepton energies within 
the whole Kf.^ Dalitz plot V and all kinematica lly allowed 
values of the Lorentz invariant x defined in ( |6.2| ). (The 
variable x determines the angle between the pion and lep- 
ton momentum for given energies Et^, Ef .) 

The inclusion of radiative corrections in this channel 
tends to give a negative shift to the phase space term. 
Assuming A4. = 0.030, and evaluating ( |7.6| ) with ao,i,2 
corresponding to a 7^ and a — Q (see Table p, one can 
check that radiative corrections induce in a negative 

shift of 1.27%. 
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7.3 Illustrative extraction of iVus 



The Kobayashi-Maskawa matrix element | Vus \ can be ex- 
tracted from the K^s decay parameters by 



\Vus\ = 



16 7r3/2 riK+ 



7r0e+z.e(7))^/' 



^+""(0)1 /(A+)i/2 
(7.22) 

The appUcation of this formula requires, of course, exper- 
imental numbers for r{K^ — > 7r'^e^i'e(7)) and A+ where 
the emission of real photons has been taken into account 
in accordance with the prescription given above. As such 
a consistent treatment of the radiative corrections cannot 
be guaranteed for all the experimental data included in 
the present values [Q of the Kes, decay parameters given 
by the Particle Data Group (PDG), the following analysis 
should only be regarded as a preliminary illustration of 
the numerics. An up-to-date determination of \Vus\ with 
the highest possible precision will then be the task of fu- 
ture high statistics measurements, and will require 
the consistent inclusion of 0{p^) terms in the form factor. 
Using § 

r{K+ -K^e+Veil)) = (2.56 ±0.03) x lO^^^MeV, 

A+ = 0.0276 ±0.0021 (7.23) 



for the parameters, S'ew(-^p, ^^z) — 1-0232 for the 
short-distance enhancement factor Q (leading logarithmic 
and QCD corrections included), and the 0{p'^) prediction 
(I7I9I) for f'^*''\Q), we find 



IK^I = 0.2173 ±0.0013 ±0.0008 ±0.0006 

= 0.2173 ±0.0016, (7.24) 

where the errors correspond to 



, , ,/ /i/+(0) 
±2-±0.05.^±^ 



iKs I ( ± 0.6% ± 0.4% ± 0.3%). 



(7.25) 



So far, our numerical estimates were based on the the- 
oretical result for the form factor up to 0{p^) in the strong 
part and to ©(e^p^) for the electromagnetic contributions. 
In the case of the pion scattering lengths, it has been 
shown explicitly that the strong interaction contribu- 
tions of 0{p^) are of comparable size to the electromag- 
netic corrections. (See also for a discussion of elec- 
tromagnetic effects in neutral pion scattering.) A similar 
feature can also be expected for the Ki^z decays. Although 
a calculation of isospin conserving corrections of 0{p^) is 
under way for the time being we have to be content 
with the rough estimate a 



given by Leutwyler and Roos already some years ago. 
Adding^ ( |7.19| ) and ( |7.26| ), we obtain 



^ (0) = 0.9874 ±0.0084, 

and, consequently, 



(7.27) 



\Vus I = 0.2207 ± 0.0013 ± 0.0008 ± 0.0019 

= 0.2207 ±0.0024, (7.28) 

which is in agreement with the current PDG value 0|, 



iKsl = 0.2196 ± 0.0023, 



(7.29) 



taken from the analysis in [^. 

On the other hand, the unitarity of the Kobayashi- 
Maskawa mixing-matrix together with the present exper- 
imental value ||l| for \yud\ implies 



IK, 



0.2287 ±0.0034, 



(7.30) 



/- 



-0.016 ±0.008 



(7.26) 



indicating a possible problem for three-generation mixing. 
Whether this discrepancy will survive once two-loop chiral 
corrections are included remains to be seen. It is interest- 
ing to notice that although these corrections would have 
to be about twice the size of the estimate (7.26) in order 
to maintain the unitarity of the CKM matrix, they still 
would not be unnaturally large from the point of view of 
the chiral expansion, and could be accounted for by 0{m?^) 
contributions Wh. 



8 Conclusions 

In the present work, we have presented general formulae 
for the Ki'i form factors at one loop in the presence of ra- 
diative corrections. Two features induced by the electro- 
magnetic interactions are worth noticing. First, the form 
factors now depend not only on the momentum transfered 
between the kaon and pion, but also on a second kinemat- 
ical variable. Second, in addition to the usual electromag- 
netic low energy constants {KiS of Ref. |0), there are 
contributions from four of the local counterterms intro- 
duced in [ p^ . These counterterms are specific to semilep- 
tonic processes of pseudoscalar mesons and renormalize 
the ultraviolet divergences induced by the exchange of vir- 
tual photons between charged mesons and leptons. 

Loops with virtual photons generate also infrared di- 
vergences. In order to deal with them, we have analysed 
the associated real photon emission processes. We have 
given a general description of the changes induced in the 
Dalitz plot density, and have proposed a model-independent 
procedure for including radiative corrections in the data 
analysis. This consists in incorporating the known long- 
distance electromagnetic effects into generalized kinemat- 
ical densities, while including all the structure dependent 
(UV sensitive) terms as corrections to the form factors. 

* Also the input parameters in the strong C'(p'*) part of the 
form factor will receive appropriate shifts j26j in a complete 
and consistent analysis of 0(j3^) effects. 
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Details of the new kinematical densities will depend even- 
tually on the way the specific experimental set-up deals 
with real photon emission. For a quite generic set-up con- 
figuration (see Sect. we have rederived explicit expres- 
sions for if mode, correcting some mistakes in earlier 
work 0. 

Within this framework, we have studied the effect of 
radiative corrections in the extraction of the CKM matrix 
element \Vus\ from the ifij mode. As compared to the 
pure 0{p'^) form factors j^], the inclusion of 0{e^p^) elec- 
tromagnetic contributions shifts /+(0) by about (0.36 ± 
0.16)%. Moreover, the radiative corrections produce an 
effective reduction of -1.27% for the phase space integral. 
We note here that the uncertainty on the form factor up 
to 0{p'^, e^p^) is well under control, and it affects the ex- 
traction of \Vus\ only marginally compared to present ex- 
perimental errors. This feature persists in the case of ifg3 
mode (not analysed numerically in this work), where the 
hadronic uncertainties are even less effective. We think 
this is a relevant result of our analysis, opening the road 
to a precision determination of |V^s|, for which the next 
two important ingredients are: 

— From the theoretical side: the inclusion of two-loop 
chiral corrections. 

— From the experimental side: a new high statistics mea- 
surement of branching ratios and slope parameters, in- 
cluding radiative corrections in the model-independent 
way outlined in this work. 
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A Mesonic Loop Functions 

The loop function HpQ{t) ||^,^ is given by 



(A.l) 



where 



/^PQ(i,M) = j^^X{t,Ml,Mi^)JpQit) 



18(47r)2~ 
-2[Ap{fi)+AQ{f,)]\, 



with 



A(a:, y, z) = + + ~ '^{xy + xz + yz), 



SpQ = Ml + MA , 



ApQ = Ml - M^, 



, , , M2 Ml 



(A.2) 

(A.3) 
(A.4) 

(A.5) 



and 



JpQ{t) 



327r2 



2 + log T7T 



Ml 



^PQ 



A 



PQ 



\^/^{t,Ml,Ml) 



X lOE 



' [t + X^/^{t,Ml,Ml)f ~ Al 



PQ 



[t- X^/^t,Ml,Ml)]^ - A 



PQ . 



1 ^Iq 



(A.6) 



Moreover, in the expansion of the form factors ,f-{t), one 
needs the function: 



KpQ{t) 



A 



PQ 



2t 



JpQ{t) 



(A.7) 



B Photonic Loop Functions 

The photonic loop contributions to the K^j, form factors 
depend on the charged lepton and meson masses m^, M^, 
as well as on the Mandelstam variables u = {pk —pif' (for 
decays) and s = (p^ +Pi)'^ (for K^^ decays). In what 
follows we denote by v the Mandelstam variable appropri- 
ate to each decay. In order to express the loop functions 
in a compact way, it is useful to define the following inter- 
mediate variables: 



R 



M2' Af2' 2Vi? 



In terms of such variables and of the dilogarithm 

dt 



Li2{x) 



t 



log(l - xt), 



(B.2) 



the functions contributing to rc{v, mj, M2; M^) and 

j.EM-loop. , . , 

jj- [V) are given by 



C{v,mj,M^ 



1 



X 



miM 1 - ^2 



-i log' X -h 2 logX log(l - ^2) - ^ -t- i log2 R 



+ Ll2{X^) + Li2{l - ^) + Li2{l - XVi?) 



(B.3) 



Riiv, mi,M^) = - - Ini? + {4 - 3Y)J^{v,mj, M'') 



r2{v,mlM^^ 



2 V u 
Ini? ^ 



T{v,mj,M^){l- R) 



{3-Y)T{v,mi,M^), (B.4) 



and 



J^{v,mlM') = -l=-^\nX. (B.5) 
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C Coefficients entering f^^'^" and /. 
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